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Abstract: We establish a general inequality on the Poisson space, yielding an upper bound 
for the distance in total variation between the law of a regular random variable with values 
in the integers and a Poisson distribution. Several applications are provided, in particular: (i) 
to deduce a set of sufficient conditions implying that a sequence of (suitably shifted) multiple 
Wiener-Ito integrals converges in distribution to a Poisson random variable, and (ii) to compute 
explicit rates of convergence for the Poisson approximation of statistics associated with geomet- 
ric random graphs with sparse connections (thus refining some findings by Lachieze-Rey and 

■ Peccati (2011)). This is the first paper studying Poisson approximations on configuration spaces 
by combining the Malliavin calculus of variations and the Chen-Stein method. 
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1 Introduction and motivation 

^ ■ The aim of this paper is to combine two powerful probabilistic techniques, namely the Chen- 

Stein method (see e.g. [21 1^) and the Malliavin calculus of variations (see e.g. [201 I25|). in 
order to study Poisson approximations for functionals of general Poisson random measures. One 

■ of our principal achievements is a general inequality on the Poisson space (see Theorem 13. ip , 
assessing the distance in total variation between the law of a Poisson random variable and the 

1/-^ \ law of a (sufficiently regular) integer-valued Poisson functional. As discussed below, a strong 

\ motivation comes from applications in stochastic geometry: in particular, our results allow to 

generalize, explain and refine some recent findings by Lachieze-Rey and Peccati [11) . dealing with 
the asymptotic behavior of general random graphs on an Euclidean space. Another remarkable 
application provided in the paper is a Poisson convergence result for sequences of 'perturbed 
multiple integrals' (see Theorem l4.1Up . extending several existing central limit theorems (CLTs) 
for multiple Wiener-Ito integrals (see e.g. [22[ l23j and [19| . respectively, for statements in the 
k> . Poisson and in the Gaussian frameworks). 

^ ■ Historically, the work [16] has been the first paper combining Stein's method for normal 

^ ■ approximations (see [J) with Malliavin calculus on a Gaussian space. This reference has been 

the seed of many generalizations and applications, for instance to fractional processes, density 
estimates and harmonic analysis of Gaussian-subordinated fields on homogenous spaces. See e.g. 
[17| for a presentation of the general theory, and [13] for several applications to the statistical 
analysis of spherical random fields. 

The present work is a natural continuation of the findings contained in [22| , where the authors 
combined Stein's method with a version of the Malliavin calculus on the Poisson space (due to 
Nualart and Vives |2U) ) in order to compute explicit bounds for CLTs involving functionals of 
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general Poissoii measures; see also |23| for analogous statements in a multi-dimensional setting. 
The two works |22| I23| have recently triggered many applications in stochastic geometry. The 
most notable papers in this respect are the following: reference j26j lays the foundations of a 
general asymptotic theory for geometric [/-statistics; reference |28| deals with CLTs for Poisson- 
Voronoi approximations; reference |29| focuses on intrinsic volumes of Poisson A;-flat processes; 
reference [7| provides an analysis of statistics associated with geometric random graphs on a 
torus; finally, reference |14j proves several multidimensional generalizations. 

In the analysis started in [26j has been extended in order to prove several necessary 
and sufficient criteria for the normal approximation of random variables having a finite chaotic 
expansion (based on the use of contraction operators - see Section BTl below). as well as to deduce 
an exhaustive asymptotic characterization of the edge-counting statistics associated with general 
stationary graphs. In the next subsection, we shall present an overview of the results established 
in that have provided the impetus for the analysis developed in the present paper. 

Remark 1.1 Other references related to the content of this paper are |10) and |18| . dealing 
respectively with the semicircular and the Marchenko-Pastur convergence in the setting of free 
probability. Note that the Marchenko-Pastur distribution is the free analogous of the Poisson 
law in classical probability theory. See also [8] for some non-standard free limit theorems, as 
well as a striking connection with the Spanish cheese industry. 

1.1 Motivation: edge counting in random geometric graphs 

Consider an integer d ^ 1, and fix the following notation: 



W :-- 



1 1 

2' 2 



d 



W:=[-l,lf, W :-- 



1 1 

"4'4 



d 



Given two mappings A — )• a(A), b{\), A > 0, we write a(A) x b(X) if there exist two positive 
constants C, C" > such that Cb{\) ^ a(A) ^ C'b{\) for A sufficiently large. For every A > 0, 
we introduce the following objects and notation: 

- rj\ IS a. Poisson measure on (Z, iF) := (W^ ^.SSiW^)), with control A£, where £ is the usual 
Lebesgue measure on (see Section [2] for precise definitions) . 

- H\ is a symmetric subset of Z x Z such that: (i) i'^{H\) < oo, and (ii) there exist H\ C Z 
such that ^ Hx and Hx = {(x, y)GZxZ: x — ye Hx}- 

- Fx = Ex,yemnw li^A y). -^A ■= and Fa := (^A - E[Fx])/ y^Y^:^. Note that the 
random variable equals the number of edges in the random graph obtained by taking 
as vertices the points of the support of r]x contained in W, and then by connecting two 
vertices x,y if and only if {x,y) G Hx- Since Hx does not contain 0, the graph obtained 
in this way has no loops, that is, no vertices of the type {x,x}. Since Hx is symmetric, 
one has necessarily that Hx = —Hx- 

- The three occupation coefficients introduced in [11, formulae (4.37)-(4.39)]: 

ip{X) := i(HxnW), Tp{X) :=£(HxnW), 4^{X) := i(HxnW). 



Remark 1.2 1. This class of geometric models contains the so-called Gilbert graphs, ob- 
tained by taking Hx = {{x,y) € M'^ x M'^ : < ||x — y\\^d < 6x}, where A — )• (5a is a 
suitable mapping with values in M+. Graphs corresponding to the cases d = 1 and d = 2 
are called, respectively, interval graphs and disk graphs- See e.g. [5l 171 [M]. 



2 



2. Since Hx can be represented in terms of the set Hx, we say that the random graph 
described above is stationary - see again 

The fohowing result is one of the main findings of [11] : it provides an exhaustive description 
of the asymptotic fluctuations of the family {-Fa}) under an additional regularity assumption on 
the occupation coefficients tpjip. 

Theorem 1.3 (See Theorems 4.9-4.11 in [111 ) Assume that x V'(A), and consider a 
standard Gaussian random variable X ^ ^(0,1). The following three properties are in order 
as A —)• oo; 

(i) If Xip{X) —?■ oo or Xip{X) x 1, then there exists a finite constant K > (independent of 
X) such that dw{Fxi^) ^ KX^'^^'^, where dwiFx^^) indicates the Wasserstein distance 
between the laws of Fx and X , so that Fx converges to X in distribution. 

(ii) If Xip{X) 0, then Var(i<\) x X'^ip{X). If in addition X'^Tp{X) oo, then there exists a 

constant K > such that dwiFx, X) ^ K ^Ay^ ip{X)^ , and consequently Fx converges 
to X in distribution. 

(iii) Assume that Xip{X) (so that Var(Fx) ^ -^^V'C'^) )> ^^'^ that the mapping X — )■ X'^ip(X) is 
bounded. If there exists a finite constant c > such that\ai[Fx) 2c, then E[F^] — > c/2, 
and F^ converges in distribution to a Poisson random variable with mean c/2 (denoted by 

Po(c/2);. 

Remark 1.4 1. For explicit examples of sets Hx verifying either one of Points (i)-(iii) in 
Theorem ll.3| see [111 Section 4.3.1]. 

2. If the constant c at Point (iii) of the previous statement is equal to zero, then one can 
prove by a direct argument that — ^ in L^. 

3. For every A > 0, the following equality in law holds: 



^A= E Wg.'A>0, (1.1) 

where 77 is a random Poisson measure with Lebesgue intensity, and Gx is a measurable 
subset of defined by the relation 

= X-^/'^Gx. (1.2) 

Point (iii) in Theorem 11.31 was proved in [11] by the method of moments, and no information 
was given about the associated rate of convergence towards the limiting Poisson distribution. 
In Section [5l we will apply the main estimates established in this paper in order to compute an 
explicit upper bound for the quantity dry (F^, Po(c/2)), where dxv indicates the total variation 
distance between the laws of F^ and Po(c/2). 

The remainder of the paper is organized as follows. Section [2] contains some preliminary 
results concerning stochastic analysis on the Poisson space. Section [3] deals with the main 
estimates proved in the paper. Section U contains an application to the Poisson approximation 
of multiple Wiener-Ito integrals. Finally, in SectionOan explicit bound for Point (iii) of Theorem 
11.31 is computed. 

Acknowledgment. I am grateful to Raphael Lachieze-Rey for many fundamental conversations 
about the topics studied in this paper. 
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2 Framework 



In what follows, the triple (Z, 2f, /i) indicates a measure space such that Z is a Borel space, 
2f is the associated Borel c-field, and /i is a non-atomic a-finite Borel measure. We define 
3ff^ := {B € Z : n{B) < oo}. Throughout the paper, we write ij = {'r]{B) : B € if^} to indicate 
a Poisson measure on (Z, 3f) with control fi. This means that r/ is a collection of random 
variables defined on some probability space indexed by the elements of and such 

that: (a) for every B,C such that B D C = 0, the random variables r]{B) and r]{C) are 

independent, and (b) for every B € i?^, r]{B) has a Poisson distribution with mean fi{B). We 
shall often write fi{B) = ri{B) — fi{B), B € if^, and fj = {fi{B) : B G The reader is referred 

e.g. to [2JJ for a general introduction to random measures of the Poisson type. By a slight abuse 
of notation, we shall sometimes write x G in order to indicate that the point x G Z is charged 
by the random measure 7/(-). 

Remark 2.1 By virtue of the specific structure of the space {Z, iF, //), we can assume through- 
out the paper that {il.,^,P) and r] are such that 



n = < 



u = ^6zj,n eNU {oo}, Zj E Z 



where 5^ denotes the Dirac mass at z, and rj is defined as the canonical mapping: 

{uj,B) ^ 7]{B){uj) = u;{B), B e 3ff,, lo e n. 

Finally, the cj-field ^ will be always supposed to be the P-completion of the c-field generated 
by r]. 

2.1 Chaos and Malliavin calculus 

Given a real p € [l,oo), the symbol U'ifJ.) is shorthand for LP(Z,2f,fi). For an integer q ^ 2, 
we shall write LP{n'^) := (Z"^ , , fi'^) , whereas stands for the subspace of L'P{n'^) 

composed of functions that are ^'^-almost everywhere symmetric. We also adopt the convention 
L'^ifJ-) = L^ifJ-) = L^ifJ-^) = -^s(m^) and use the following notation: for every q ^ 1 and every 

1/2 

J K^l, ■■■,^q}yK-^l, ^qJH'-KU^l, ■■■,"'^qj, \\J ||L^(^j9) — \J , J / 

Z1 



{f,9)L2(t^l) = / fizi,...,Zq)g{zi,...,Zq)n'l{dzl,...,dZg), ||/||l2(|.'J) = (/,/)/2(^,)- 



For every / € L^(/u'^), we denote by / the canonical symmetrization of /. Plainly, ||/||i2(-^tj') ^ 

ll/llL2{/i9). 

Definition 2.2 For every deterministic function h G L'^{n), we write 

h{h)=fj{h) = [ h{z)f,{dz) 
JZ 

to indicate the Wiener-Ito integral of h with respect to f). For every q ^ 2 and every / € L^(/u'^), 
we denote by Iq{f) the multiple Wiener-Ito integral, of order q, of / with respect to f/. We also 
set Iq{f) = Iqif), for every / G L^(//'') (not necessarily symmetric), and lo(^) = b for every real 
constant b. 
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The reader is referred for instance to |211 Chapter 5] for an exhaustive discussion of mul- 
tiple Wiener- Ito integrals and their properties (including the forthcoming Proposition 12.31 and 
Proposition I2.4|) . 

Proposition 2.3 The following equalities hold for every q,m^ 1, every f G L1{^'^) and every 
1. E[I,{f)] = 0, 

2- E[Iq{f)Im{g)] = ql{f,g)L^(^i)l(q=m) (isomctric property). 

The Hilbert space composed of the random variables with the form Iq{f), where q ^ 1 and 
/ G (//''), is called the gth Wiener chaos associated with the Poisson measure r]. The fol- 
lowing well-known chaotic representation property is an fundamental feature of Poisson random 
measures. Recall that ^ is assumed to be generated by rj. 

Proposition 2.4 (Wiener-Ito chaotic decomposition) Every random variable 

F G L^{n,^,P) := L^{P) 

admits a (unique) chaotic decomposition of the type 
oo 

F = E[F] + Y,im, (2.3) 

i=l 

where the series converges in L^{P) and, for each i^ 1, the kernel fi is an element of L1{^'^). 

For the rest of the paper, we shall use definitions and results associated with the Malliavin- 
type operators defined on the space of functionals of the Poisson measure 7/. Our formalism 
coincides with the one introduced by Nualart and Vives in [20] • A more recent introduction 
to the Malliavin calculus on the Poisson space can be found in Privault's monograph |25| : for 
the convenience of the reader, some basic definitions and results are presented in the Appendix 
(see Section [6|). In particular, we shall denote by Z), 5, L and , respectively, the Malliavin 
derivative^ the divergence operator, the Ornstein-Uhlenbeck generator and its pseudo-inverse. 
The domains oi D, 6 and L are written domD, douiS and domL. The domain of L^^ is 
given by the subclass of L'^{P) composed of centered random variables. Given a not necessarily 
centered F G -L^(P), one sets by convention L'^F = L-^{F-E{F)), so that LL'^E = F-E{F) 
for every F G L'^{P). Since the underlying probability space is assumed to be the collection 
of discrete measures described in Remark 12.11 one can meaningfully define the random variable 
u I— 7- Fz{ijj) = F{ljj + 5z), w G Jl, for every given random variable F and every z Z , where 6z 
is the Dirac mass at z. One therefore has the following neat representation of D as a difference 
operator: 

Lemma 2.5 For each F G dovaD, 

DzF{io) = Fz{uj) - F{uj), a.s.-P{dio), a.e.-fi{dz). 
A proof of Lemma 12.51 is given in [20j . 
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3 A general inequality on the Poisson space 



The following result is the main finding of the paper. As anticipated in the Introduction, 
the proof makes use of the so-called Chen-Stein method for Poisson approximations. A classic 
reference on the subject is the book by Barbour et al. [2]; more recent references are the two 
surveys [SI'S]. Recall that the total variation distance between the laws of two random variables 
X, Y with values in := {0, 1, 2, ...} is given by 



dTv{X,Y)= sup \P{X eA)-P{Y &A)\ = ]^^\P{X = k)-P{Y = k)\. (3.4) 

f is strictly 



Of course, the topology induced by dxv on the class of all probability laws on 
stronger than the topology induced by convergence in distribution. 

Given a function / : — )> M, we denote by A/ the forward difference given by Af(k) := 
f{k + 1) — f{k), k = 0,1,2,...; we also use the symbol A^/ = A(A/). Finally, we write 

ll/lloo = SUPfcg2;+ 1/(^)1- 

Theorem 3.1 (Malliavin bounds for Poisson approximations) Fix c > 0, and let Po(c) 
indicate a Poisson random variable with mean c. Assume that F G L'^{P) is an element of 
domZ? such that E{F) = c and F takes values in Z+. Then, 



dTv{F,Po{c)) 



1 



-E\c-{DF,-DL-'F)l2^^) 



-E 



\D,FiD,F-l)D,L-^F\fiidz] 



E 



+- 



-E 



-{DF,-DL-^F)L2^^)y 
ID^F (D^F - 1) D^L-^F\ n{dz] 



(3.5) 
(3.6) 

(3.7) 



Remark 3.2 1. Let i*" be a centered element of domi^ such that E{F'^) = 1, and let Z ^ 
^(0, 1) be a standard Gaussian random variable. In |22| Theorem 3.1], it is proved that 



dwiF,Z) ^ E\l- {DF,-DL-'F)l2 



''Ml 



+E 



I \{D,FfD,L-^F\fi{dz) 

JZ 



(3.8) 
(3.9) 



where d^y denotes the Wasserstein distance between the laws of two random variables. 

2. If F = 7?(^), where /i(^) = c, then one has that D^F = —DzL^^F = and therefore 

c-{DF,-DL-^F)l2u)= I \D,F{D,F-l)D,L-^F\f,{dz) = 0, 

Jz 

as expected. 

3. Let c, c' > 0. Standard computations (see e.g. [U Corollary 3.1]) yield that 

dTviPo{c),Fo{c)) i:\c- c\. (3.10) 

It follows that the content of Theorem 13.11 can be extended to a random variable F with 
arbitrary positive expectation E{F) = d by using the triangular inequality: 

dTy(F,Po(c)) ^dTv(Po(c'),Po(c)) + dTy(F,Po(c')) ^ \c-d\+dTv{F,Vo{c')). 
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4. Integrating by parts, one sees that E[{DF, —DL F)^2(^)] = Var(F). 

Proof of Theorem For every A C Z+ , we denote by : Z+ R the unique solution to the 
Chen-Stein equation 

lA(k) - P(Po(c) £A) = cf{k + 1) - kf{k), A: = 0, 1, 

verifying the boundary condition A^/(0) = 0. Combining e.g. [91 Theorem 2.3] with |6i Theorem 
1.3 and pp. 583-584], we immediately deduce the estimates 

f^lloo ^min( , ||A/a||oo ^ ^ and ||AVa||oo ^ ^ "J"" • (3.11) 



Using the relation 5D = —L one infers that 

P(Po(c) eA)- P{F e A) = E[FfA{F) - cfA{F + 1)] 

= E[{F - c)fA{F) - cAfAiF)] = E[6{-DL-'F)fA{F) - cA/a(F)]. 



Integrating by parts yields 

E[6{-DL-'F)fA{F)] = E[{DfA{F),-DL-'F)L2 



where, by virtue of Lemma 12.51 DzfA{F) = fA{F + D^F) — fA{F). Observe that Lemma 12.51 
implies that, since F takes values in Z_|_, then one can always choose a version of D^F with 
values in Z, in such a way that F + D^F = F^ takes values in Z_|_. Now, for every / : Z_(_ — >■ R 
and every k,a £ TLj^ such that k > a, one has that 

fc-i 

fik) = f{a) + A/(a)(fc -a)+Y, A^f{j){k - 1 - j); 

j=a 

on the other hand, when k,a £ Z4. are such that k < a, 

a-l 

fik) = f{a) + Af{a){k -a) + Y, ^^f{j){j + l-k). 

j=k 

These two relations yield that, for every k,a £ Z_|_, 

IIA^ f II 

\f{k) - f{a) - Af{a){k - a)\ ^ \{k - a){k - a - 1)\ . 

Taking a = F and k = F^, one therefore deduces that 

DjA{F) = AfA{F)D,F + R,, 
where Rz is a residual random function verifying 

\Rz\^^^^^\DzF{DzF-l)\, zeZ. 

As a consequence, 

P(Po(c) eA)- F{F £A) = E[{DfA{F), -DL-'F)l2^^) - cAfA{F)] 

= E[AfA{F){{DF, -DL-^F)l2^^) - c)] - E f {R, x D,L-^F) i^{dz), 

Jz 

and the desired conclusion follows by taking absolute values on both sides, as well as by applying 
the estimates (|3.1ip . Inequality (|3.7p follows the Cauchy-Schwarz inequality. 
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□ 

The following statement is an immediate consequence of Theorem 13.11 

Proposition 3.3 (Poisson limit theorems) Fix c > 0. Let {F„ : n ^ 1} C AovnD be a 
sequence of random variables with values in Z_|_ such that — t- c, as n oo. Assume that, 
as n oo, 

1. E\c- {DFn,-DL~^Fn)L^^,)\ ^0, and 

2. E[J^\D,FniD,Fn-l)D,L-^Fn\fJ.idz)] ^0. 

Then, limn_!.oo c^Ty(-^n) Po(c)) — )• and F^ converges in distribution to Po(c). 

Remark 3.4 It is well known that the Poisson distribution is determined by its moments (see 
e.g. \2ll p. 42-43]). It follows that, in order to prove that a given sequence {Fn} converges 
in distribution to Po(c), it is sufficient to prove that E[F^] — )• £'[Po(c)'^], for every integer 
k ^ 1. This is the so called 'method of moments', requiring a determination of all moments 
associated with each Fn- Although popular in a Gaussian setting (see [21] for an overview), 
such a technique is extremely demanding (and very little used) in the framework of Poisson 
measures. This is mainly due to the fact that the combinatorial structures involved in the so- 
called 'diagram formulae' (that are mnemonic devices used to compute moments by means of 
combinatorial enumerations - see [21', Chapter 4]) become quickly too complex to be effectively 
put into use. One should compare this situation with the statement of Proposition 13.31 which 
only involves two sequences of mathematical expectations. In the forthcoming Sections [JHSl we 
will see that, in many applications, in order to check Points 1 and 2 in Proposition 13.31 one is 
naturally led to assess and control expressions that are only related to the second, third and 
fourth moments of the sequence {Fn}. 

4 Applications to multiple Wiener-Ito integrals 

In this section, we use Theorem 13. II in order to establish Poisson convergence results for general 
sequences of Z+-valued random variables of the type 

Fn = Xn + Bn + Iq{fn), n ^ 1, (4.12) 

where: (i) {x„ : n ^ 1} is a sequence of positive real numbers, (ii) (7 ^ 2 is an integer independent 
of n, (iii) Iq indicates a multiple Wiener-Ito integral of order q, with respect to the compensated 
measure fj, where r/ is a Poisson measure on the Borel space (Z, 3f) with cr-finite and non-atomic 
control measure /i, (iv) /„ G L'^{fi'i), and (v) {Bn : n ^ 1} is a smooth vanishing perturbation, 
in the sense of the forthcoming Definition 14.21 

Remark 4.1 1. In the statement of the forthcoming Theorem 14. 101 we shall implicitly allow 
that the underlying Poisson measure r] also changes with n. In particular, one can assume 
that the associated control measure fi = fin explicitly depends on n. As discussed in 
Section [5l this general framework is useful for geometric applications. 

2. We consider perturbed sequences of multiple integrals because, in general, it is not clear 
whether a non-trivial random variable of the type x + Iq{f), where x € M_|_ and q ^ 2, can 
take values in Z4., and therefore whether Theorem 13.11 can be directly applied. However, 
as seen e.g. in Section [5l in applications one often encounters sequences of integer- valued 
random variables whose chaotic decomposition is such that all terms except one vanish 
asymptotically. 
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Definition 4.2 (Smooth vanishing perturbations) A sequence {Bn : n ^ 1} C L'^(P) is 
called a smooth vanishing perturbation if Bn, L~^Bn G domD for every n ^ 1, and the following 
properties hold: 



lim E[BI] = 



lim E 

n— >cxD 

lim E 



\DBr, 



hm 
hm E 

n— )-oo 



0, 
0. 



(4.13) 
(4.14) 
(4.15) 



Note that, if (|4.14p - (|4.15p are verified, an application of the Cauchy-Schwarz inequality yields 
that 



lim E 

n— )-oo 



\DB. 



n|lL3(/i) 



lim E 

n— >oo 







Remark 4.3 Using a Mehler-type representation of the Ornstein-Uhlenbeck semigroup (such as 
the one stated in |[25i Lemma 6.8.1]), one sees that the following inequalities are always verified: 



E 



\DB.^ 



^ E 



E 



\DB,. 



^ E 



Remark 4.4 The following conditions (a)-(b) are sufficient in order for a given sequence {Bn} 
to be a smooth vanishing perturbation. 

(a) There exists an integer M ^ 1, independent of n, such that 



M 



M 



Bn = Ii{gi,n) ■■= Bi^n, n ^ 1, 



i=l 



i=l 



where € L^(//*) for every i = 1,...,M. Note that these assumptions imply that 
E[Bn] = and Bn,Bi^n € domD for every n and every i. 



(b) As n — 7> oo, for every i = I, M, 



E[Bl]^0, E 



{D.B.^nffiidz] 



0, and E 



[ {D,B,^nfii{dz) ^ 0. (4.16) 
Jz 



Explicit examples of smooth vanishing perturbations appear in Section [5] - see in particular 
Remark 15.21 



4.1 Digression: stars and products 

In order to state and prove the main results of this section, it is now necessary to introduce 
contraction operators and to discuss their role in product formulae. 

The kernel f -k^ g on z^^'^^'^^\ associated with functions / G L^^fi^) and g € L'^{fi'^), where 
p,q ^ 1, r = 1, . . . ,p A q and / = 1, . . . , r, is defined as follows: 

, ^r—h ^1) ) • • • ) tp—n Si) ) • • • ) Sq—r) (4.17) 
,dzi)f{zi, , . . . , Zi,^i, . . . ,7r_/,tl, , . . . ,tp-r) 

Xg{zi,, ... ,2;/, 71, . . .,Jr-l,Sl, ,Sq-r). 



H^dzi, ... 
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As it is evident, the operator ' ' reduces the number of variables in the tensor product of / and 
g from p + q to p + q — r — I: this reduction is obtained by first identifying r variables in / and 
g, and then by integrating out / among them. To deal with the case / = for r = 0, . . . ,p A q, 
we set 

/ *r- dilli ■ ■ ■ 1 ^r-i tl, 1 ■ ■ ■ 1 tp—r: Si, , . . . , Sq—r) 
= fin, • • • ,7r,*l, , • • • ,tp-r)g{7l,- ■■,7r,Sl,,... ,Sq-r), 

and 

/*0 9{tl, ,...,tp,Sl,,...,Sq) = f (^g{ti,,.. .,tp,Sl, ,...,Sq)= f{ti, , . . .,tp)g{si, , . . .,Sq). 

By using the Cauchy-Schwarz inequality, one sees immediately that f g is square-integrable 
for any choice of r = 0, . . . ,p A g , and every / G L^fi^), g G L^(/i''). 



Remark 4.5 For every 1 ^ p ^ q and every r = 1, ...,p, 

if 4 gfdf.P+'i-^ = [ if f)ig g)di/, (4.18) 



L 



ZP+q-r 

for every / G L^fiP) and every g G L'^i^fJ-' 

The next result is a fundamental product formula for Poisson multiple integrals (see e.g. 
for a proof). 

Proposition 4.6 (Product formula) Let f G L1{ijP) and g G L^fi'^), p,q^ I, and suppose 
moreover that f -k\. g ^ L^(/iP+'^^''^') for every r = 1, . . . ,p A q and I = 1, . . . ,r such that I ^ r. 
Then, 

W)hi9) = E ( f ) ( ' ) E ( l) Wr-^ (M^) , (4.19) 

with the tilde ~ indicating a symmetrization. 

In order to be able to directly apply the computations contained in |22l Proof of Theorem 
4.2], in what follows we shall always work under the following technical assumption. 

Assumption 4.7 (Assumption on integrands) Every random variable of the type Y = 
Iq{f), where q ^ 2 and / G L'^n'^), considered in the sequel of this paper is such that the 
following properties (i)-(iii) are verified. 

(i) For every r = 1, ...q, the kernel /j *q~^ fi is an element of L'^{i/). 

(ii) Every contraction of the type [zi, ...,Z2q~r-l) ^ *r \fi\{zi, Z2q~r-l) is well-defined 
and finite for every r = 1, g, every / = 1, r and every (^i, Z2q-r-l) ^ 

(iii) For every /c = 1, 2q — 2 and every (r, I) verifying k = 2q — 2 — r — I, 
{f{z,-)-k\. f{z,-)Y dn^ fi{dz) <Qo, 

where, for every fixed z G Z, the symbol f{z,-) denotes the mapping (zi, Zq-i) i-^ 

f{z,Zl, ...,2g-l). 



10 



Remark 4.8 According to the discussion contained in [22], Point (i) in Assumption 14.71 implies 
that the following properties (a)-(c) are verified: 

(a) for every r = 1, ...,q and every / = 1, r, the contraction / ★J. / is a well-defined element 
of 

(b) for every r = 1, q, f *r f is an element of L^(^^'^~'^'); 

(c) for every r = 1, ...,q, and every / = 1, ...,r A {q — 1), the kernel / / is a well-defined 
element of L^(//^'^~''~'). 

In particular, the multiplication formula ()4.19p implies that every random variable Y veri- 
fying Assumption 14.71 is such that Y'^ G L'^{P), yielding in turn that £'[y^] < oo. Analogously, 
one can also show that, under Assumption 14. 7| the random variable {DY, — D L^'^Y) ^2 ^^^-^ is 
square-integrable (and not merely an element of L^{P)). 

Remark 4.9 For instance. Assumption 14.71 is verified whenever / is a bounded function with 
support in a rectangle of the type B x ■ ■ ■ x B, where fJ.{B) < oo. 

4.2 Poisson limit theorems 

The following statement contains the main result of the section. 

Theorem 4.10 (Poisson limit theorems on a perturbed chaos) Fix c > and let Po(c) 
be a Poisson random variable with mean c. For a fixed q ^ 2, let {F„ : n ^ 1} be a "Lj^-valued 
sequence as in ^.12^ , such that Xn c and E[Iq{fn)^] — )• c. Assume moreover that the following 
three conditions hold: 

(i) For every n ^ 1, the kernel fn verifies Assumption \4-'i\ 

(ii) For every r = 1, ...,q, and every I = I, ...,rA{q — l), \\ fn*i fn\\L'^(p2q-r-i^ ^ (as n ^ oo). 

(iii) The relation sup„ ||/n||L4(^q-) < oo holds and, as n ^ oo, 

I {f^ + q\'f^-2q\f^)dfi''^0. (4.20) 
Jzi 

Then, (iry(F„, Po(c)) — > 0, as n oo. 

Remark 4.11 1. Condition (|4.20p is trivially verified whenever 

fn{zi, ■..,Zq) = ^ iHnizi, ■■■,Zq), 

where Hn is some measurable symmetric subset of Z'^. 

2. When q = 2, Conditions (ii) and (iii) in the statement of Theorem 14 . 1 1 boil down to the 
following asymptotic relations 

ll/n *1 fn\\L^{f,2) 0, Wfn +2 fuWl^i^^) 0, and 
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Theorem 14. 101 should be compared with the following central limit result, first proved in [221 
Theorem 4.2]. 

Theorem 4.12 (CLTs on a fixed chaos, see [22]) Let X ~ ^(0,1). Fix q ^ 2, and let 

Fn = Iq{fn), n ^ 1, be a sequence of multiple stochastic Wiener-Ito integrals of order q. Suppose 
that, as n ^ oo, E{F^) — ?> 1. Assume in addition that the following conditions are verified: 

(i) For every n ^ 1, the kernel fn verifies Assumption \4-7\ 

(ii) For every r = 1, ...,q, and every I = 1, ...,rA{q — l), ||/n*r /n|lL2(/i2<7-r-i) (^as n —>■ ooj. 

(iii) As n ^ oo, 

[ f^dfi'^^0. (4.21) 
Jzi 

Then, Fn converges in distribution to X, as n ^ oo, in the sense of the Wasserstein distance. 



Proof of Theorem \4-l(J\ We have to prove that, if (i)-(iii) are verified, then Point 1 and Point 2 
in the statement of Proposition 13.31 hold. The proof is divided into three steps. 

Step 1: Computations related to DIq{fn)- One has that Dzlq{fn) = Q^g-iifniz, ■) and 

{DIqifn),-DL-'lqifn))LH>.) = ^\\DIqifn)\\LH>.)- 

Using the computations contained in [22', p. 464], one sees that 

2q-2 

{DJqifnf} = q'Yl IpiGl-'fiz, •)), (4.22) 
p=0 

where 

Gf'f{z,-){zi,...,Zp) = ^^1^2,-2-r-i=p}r\i\ I (M/(z,.)4/(^r)(^i,-,^p), (4.23) 

and the stochastic integrals are set equal to zero on the exceptional set composed of those z such 
that f{z, •) f{z, •) is not an element of L^(/i^''"^^''~') for some r, /. Combining the triangular 
and Cauchy-Schwarz inequalities with the estimates in \22\ Theorem 4.2] (see also [111 Theorem 
3.5] for a more compact statement), one deduces that there exists a constant K, depending 
uniquely on c and q, such that, writing E[Iq{fn)'^] ■= Vn, 

E[\c-q-^\\DIq{fn)\\l2^^)\]^\c-yn\ + K X max ^ Wfn 4 fn\\L^^.^^-r-') (4.24) 

l=l,...,rA{q-l) 

^E[\\DIq{fn)\\\,^^^] ^ K I ^max^ + ||/n|li4(^) \ . (4.25) 

b=l,...,rA('<?-l) J 

In particular, these relations imply that the sequence n i-)- E[\\DIq{fn)\\^^k(^^-^] is bounded for 
k = 2,3,4. 
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Step 2: Dealing with Bn- Since {-Bn} is a smooth vanishing perturbation, using the Cauchy- 
Schwarz inequahty and the results from Step 1, one deduces immediately that the sequence 

{DFn,-DL-'Fn)L^^,) - {DI,{U),-DL-'l,{fn))LH^) 

= (Dlgifn), -DL-^Bn)L^f.) + {DBn, -DL-^Ig{U))L^^.) + {DBn, -DL-^Br,)L2^^) 

converges to zero in L^iP). Exploiting in a similar way the fact that {Bn} is a smooth vanishing 
perturbation together with the bounded character of |1|2(^)] and E[\\DIq{fn)\\^4^^^], 

one also deduces that 



E 



[ \D,Fn {D,Fn - 1) D,L-^Fn\ li{dz) 
Jz 

- [ \DJg{fn) (DJgifn) - 1) D,L-^Ig{U)\fl{dZ 

Jz 

Now observe that, thanks to (|4.24|) and Assumption (ii) in the statement, 

E[\c - q-'\\DlM\\l2^^)\] ^0, n ^ oo. 
It follows that the proof is concluded once we prove that, as n ^ oo. 



0. 



E 



[ \DJg{fn) {DJgifn) - 1) D,L-^Ig{fn)\ 

Jz 



q-^E 



0. 



This is the object of the forthcoming Step 3. 

Step 3: Proving \4.26 ). Using the Cauchy-Schwarz inequality, we infer that 



E 



[ {DJg{fn)f \DJgiU) - iHdz) 

JZ 



(4.26) 



^ ^E\\Dlg{U)\\\.(^^^ V X ^ P./g(/n) - 1)2] ll{dz) 

Also, recall that (^G\ ^ fn{z, •) = qq\ f^g-i ■)'^ dfi'^^^ , and consequently 

i^[||^^a(/n)lli2(^)]=gg!||/n|li2(^,). 

From (|4.22p . we infer that 

E [{DJg{fn)f {DJgUn) - If] = (1 + q^Gl^Uz, •)) X q^Gl^Uz, •) 

+q' p\\\Gr'f{z,-))\\h^,.) 

p=l,...,2q-2 

+q% -ly.j^ GlZ\Uz, ■) (gI^IUz, •) - 2q-'fiz, •)) 

Integrating and simplifying the RHS of the previous equality by reasoning as in |22^ p. 467], one 
deduces from Assumption (ii) in the statement that 

E [{DJg{fn)f {DJgiU) - If] is{dz) = oil) + / {qq\f^ + qq\^f^ - 2qq\^f^) d/x"?, 

Jzi 

where o(l) indicates a sequence converging to 0, as n ^ oo. The conclusion follows by applying 
relation (g^O]). □ 
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5 An application to geometric random graphs 

The following statement provides the announced refinement of Theorem 11.31 

Theorem 5.1 Let the notation and assumptions of Theorem M . 'J[ - ( Hi) prevail. Then, there exists 
a finite constant K , independent of X, such that 



dTviF^, Po(c/2)) ^ \E[F^] - c/2\ + K,/Mf{X). (5.27) 
Proof. Using e.g. [261 Theorem 3.1], one sees that admits the following chaotic decomposition 

F^ = E[F^]+h{fi,x) + l2{f2,x), 

where Ii , I2 indicate (multiple) Wiener-Ito integrals with respect to the compensated measure 
Vx = V- fi,x{z) = >^fw ^H^{z,x)£{dx)lw{z), and f2,x{zi,Z2) = ^'i-Hxn{WxW){zi, Z2)- Also, 
Campbell's Theorem (see [27, Theorem 3.1.3]) implies that = ^X'^£{Hx). Applying The- 

orem [3TT] together with Remark 13.21 and the Cauchy-Schwarz inequality yields that 

dTviFlPoic/2)) < \E[F^^] - c/2\ + ^ ~ So(A) + i^ijl!^Hi(A) x E^W, 

where 



)(A) 



KA) 



E 



{E[F*] - {DF;^,-DL-^F*)L2^,,)y 



Si(A) := JXE 



{D.L-^F^ 



1 P*\2l 



XE \ [ {D,F*)'^{D,F* - lyeidz) 



One has that D^F^ = fi^x{z) + 2Ii{f2^x{z, •))• Using the computations contained in |22| p. 470], 
we infer that 

{DFl-DL-'Ft^)L2^x,^ 

= II/i,a|Il2(A£) + 2|I/2,a|Il2({A^)2) + 2A(/2,A *2 /2,a) + 2/2(/2,A *1 f2,x) + 3/l(/l,A *2 f2,x)- 

Using the computations contained in [lT| Theorem 4.8] one sees that the following five relations 
are verified: 

(1) \\fi,x\\h^,,)^X^i,{Xf^Xi,{X), 

(2) 2||/2,A|li2((A£)2) = E[F^], 

(3) |I/2,A*^ f2,x\\h^x,^ ^ A3^(A)2 X A^(A), 

(4) ||/2,A*{ /2,A|li2((,,)2) ^ AV(A)3 X (A^(A))2, 

(5) ll/i,A*{ f2,x\\h^xe) - A5^(A)4 X (AV'(A))3. 
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This implies in particular that Ho(A) x \J Ai/;(A). By using the explicit expression —D^L^^F^ = 
fi,x{z) + -^i(/2,a(-2, •))' it is not difficult to prove that the mapping A i— t- Hi(A) is necessarily 
bounded, so that the statement is proved once we show that H2(A) x y^AV'(A). U sing the 
relation 

(D.F^f = \fi,xizf+iX [ f2,xiz,x)iidx) 

JZ 

+4/i[/2,a(^, ■){hM + /2,a(^, •))] + 4/2(/2,A A /2,a), 

developing the square {DzF^ — 1)^ and integrating with respect to z yields 



XE 



[ {D^FtfiD, 

JZ 



n - ^f^idz) 



A / ifiAz)' + h,x{^)yidz)+8X' 




z Jz 



+48A^ 



+16A' 



z \Jz 



fl,xiz)f2,xiz,xflidx)i{dz) 

f2,x{^,x)-'£{dx)] e{dz)+AX' I [ h,x{z,xfe{dx)e{dz) 



z Jz 



Z JZ 



flxi^, X){h,x{z) + /2,a(^, ^))(/i,a(^) + /2,a(^, x) - mdx)l{dz). 



Exploiting the explicit form of /2,A) one sees that 
4A^ / f2,x{z,xfl{dx)l{dz) + lQX^ [ f2,xiz,x)^iidx)i{dz)-16X'^ [ f2.x{z,xfl{dx)l{dz)={). 

Jz ' JZ ' JZ ' 

One can now deal with the remaining terms by using once again the estimates contained in 
[TTl Theorem 4.8]: this entails the relation XE[j^{DzF^f{D^F^-lfl{dz)\ x XiIj{X), and 

□ 



consequently the desired conclusion H2(A) x ^JX'^l){X). 



Remark 5.2 By inspection of the previous proof, one sees that the mapping A i-^ A(/i,a) is a 
smooth vanishing perturbation, in the sense of Definition 14.21 It follows that the convergence 
Po(c/2) can alternatively be proved by directly applying Theorem 14.101 
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6 Appendix: Malliavin operators on the Poisson space 

We now define some Malliavin-type operators associated with a Poisson measure r], on the Borel 
space {Z,2f), with non-atomic control measure ^. We follow the work by Nualart and Vives 

m- 

The derivative operator D. 

For every F € L'^{P), the derivative of F, DF is defined as an element of L^{P; L'^{h)), 
that is, of the space of the jointly measurable random functions u : il. x Z M such that 
E [Jr^ulfiidz]] < oo. 

Definition 6.1 1. The domain of the derivative operator D, written domD, is the set of all 
random variables F G L'^{P) admitting a chaotic decomposition (OP such that 

^kk\\\fk\\\2(^^k) < oo, 
2. For any F G domD, the random function z D^F is defined by 

oo 

D,F = ^kh-i{fk{zr))- 

The divergence operator 6. 

Thanks to the chaotic representation property of ry, every random function u € L'^{P, L^(/u)) 
admits a unique representation of the type 

oo 

Uz = Y,Wk{z,-)), zeZ, (6.28) 

where the kernel is a function of A; -|- 1 variables, and fk(z, •) is an element of Ll{ix^). The 
divergence operator 6{u) maps a random function u in its domain to an element of L'^{P). 

Definition 6.2 1. The domain of the divergence operator, denoted by dom5, is the collection 
of all u € {P, L'^ {n)) having the above chaotic expansion 46'. 28\) satisfied the condition: 

k^O 
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2. For u S doraS, the random variable 5{u) is given by 

where f}^ is the canonical symmetrization of the k +1 variables function f^. 

As made clear in the following statement, the operator 5 is indeed the adjoint operator of D. 
Lemma 6.3 (Integration by parts) For every G € doraD and u G doxaS, one has that 

E[G5{u)]=E[{DG,u)l^^,)\. 
The proof of Lemma [6. 31 is detailed e.g. in |20] . 

The Ornstein-Uhlenbeck generator L. 

Definition 6.4 1. The domain of the Ornstein-Uhlenbeck generator, denoted by domL, is 
the collection of all F £ L^{P) whose chaotic representation verifies the condition: 

2. The Ornstein- Uhlenbeck generator L acts on random variable F G domL as follows: 

LF = -J^klkifk). 

The pseudo- inverse of L. 

Definition 6.5 1. The domain of the pseudo-inverse of the Ornstein-Uhlenbeck generator, 
denoted by L^'^ , is the space Lq(P) of centered random variables in L^{P). 

2. ForF=j: hifk) G Ll{P) , we set 

L-'F = -Y.\h{fk). 
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